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INTRODUCTION 


'I'he  spt-v  iral  raiiiuiue  of  the  sky  must  be  determiiuxl  in  order  to  predict 
ttu‘  perforiiuince  of  passive  infrared  surveillance  systems  operating  against  air¬ 
borne  targets  Spectral  radiance  is  esp<<ially  important  when  viewing  targets 
close  to  the  horizon  For  these  low  angles  of  view  ing,  the  smooth  sea  is  a  poor 
emitter  of  infriired  wavelengths,  and  its  radiance  is  primarily  determined  by  the 
eniis-sions  of  the  intervening  atmosphere  and  reflections  from  the  sky  For  a 
w  ind-rutllwi  sea,  the  eia.ssions  from  the  individual  wave  facets  must  also  be  ac¬ 
counted  for  in  the  background  scene  In  the  spectral  regions  about  the  6  3  ^m 
water  vapor  band,  where  the  aljsoriition,  and  consequently  the  emissivity,  is 
high,  the  sky  is  (*ss*'ntially  a  blacklKaiy  at  the  temperature  of  the  lowest  atmos¬ 
pheric  layer,  with  little  variation  with  elevation  angle  In  the  8-12-|im  region, 
the  alisorption  (and  emissivity)  depends  on  the  optical  path  length,  and  therefore 
the  radiance  will  1m-  strongly  dependent  on  the  elevation  angle  (!loud  ty]H'  and 
coverage  w  ill  also  have  strong  efr**cLs  on  radiance,  depe  nding  ujMin  the  elevations 
and  temperatures  of  the  clouds 

In  the  absence  of  radiometric  measurements,  we  must  pr»*sently  rely  on 
the  LOWTllAN  atmospheric  transmittance/radiance  computer  code  (Kneizys,  et 
al  .  lUiS3),  along  with  measured  meteorological  parameters  and  models  of  aerosol 
size  distributions,  to  predict  the  infrared  background  emi.ssions  The  Lf)WTRAN 
6  code  hits  proven  to  be  a  versatile  tool  in  predicting  atmospheric  radiance  above 
the  horizon  (Hughes  and  Jensen,  1988).  A  model  has  recently  been  developed  of 
the  effective  radiance  of  the  sea  surface  as  a  function  of  the  viewing  zenith  angle 
and  sensor  height  fWollenweber.  1988b).  and  has  been  incorporated  into 
LOWTRAN  6  f^or  processing  on  an  HP-9020  computer.  The  model  (see  Appendix 
A)  uses  the  wave  slope  statistics  of  Cox  and  Munk,  1954,  which  is  based  on  sur¬ 
face  wind  speetJ  and  direction  relative  to  the  look  angle.  The  model  is  currently 
limited  to  the  Navy  Maritime  Aerosol  Model  (see  Appendix  B),  which  in  addition 
to  the  standard  meteorological  parameters,  requires  as  inputs  the  current  and 
24-hour  averaged  wind  speeds  and  the  origin  of  the  air  mass.  Basically,  the 
NOSC-developed  model  calculates  the  total  contributions  to  radiance  at  the  sen¬ 
sor  from  the  emissions  of  the  intervening  atmosphere  and  the  sky  reflections  and 
surface  emissions  from  the  wave  slope  surfaces.  To-date  this  model  has  received 
very  little  validation,  and  the  computer  code  has  not  yet  been  made  available  to 
the  scientific  community. 

In  this  report,  measurements  of  near-horizon  infrared  (8-12-|im)  radiances 
made  with  a  calibrated  thermal  imaging  system  are  presented.  The  radiance 
“scenes”  (thermograms)  are  used  to  compare  the  mean  equivalent  blackbody  tem¬ 
peratures  corresponding  to  an  area  1  degree  above  and  1  degree  below  the  hori¬ 
zon  during  different  meteorological  and  surface-wind-speed  conditions.  Airborne 
measurements  of  the  vertical  profiles  of  meteorological  parsuneters  and  sea  sur¬ 
face  temperatures  £ire  then  used  to  evaluate  the  ability  of  the  radiance  algorithms 
with  the  LOWTHAN  6  code  to  predict  the  measured  results. 

MEASUREMENTS 

For  this  study,  a  Piper  Naveyo  aircraft,  equipped  with  Rosemount  tem¬ 
perature  and  pressure  probes  and  an  EG&G  dewpoint  sensor,  made  vertical 
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spirals  over  the  oceein  to  obtain  temperature,  relative  humidity,  and  pressure, 
which  are  required  inputs  to  the  LOWTRAN  6  computer  code  for  calculating  the 
sea  and  sky  radiances.  A  Bsu’nes  PRT-5  radiation  thermometer  was  also  onboard 
the  aircraft  to  measure  the  actual  sea  surface  temperatures  from  low-level 
( <  20m),  constant-altitude  flights.  At  the  time  the  meteorological  parameters 
were  obtained,  measurements  of  infrared  (8-12-nm)  horizon  radi6mces  were  also 
made  with  a  calibrated  thermal  imaging  system  (AGA  THERMOVISION,  Model 
780)  and  a  2.95-degree  field-of-view  lens. 

For  these  measurements,  the  scanner  was  located  at  an  elevation  of  33  m 
on  the  Point  Loma  peninsula  in  San  Diego  and  was  directed  westward  over  the 
ocean.  The  response  of  the  system  is  determined  by  placing  a  blackbody  of 
known  temperature  (±0.1  °C  for  temperatures  <50°C)  in  front  of  the  lens  aper¬ 
ture.  The  digitized  video  signal  transfer  function  of  the  system  then  allows  the 
blackbody  temperature  to  be  reproduced  to  within  ±0.2  °C. 

Examples  of  measured  radiance  scenes  are  shown  in  Fig.  1  and  2  for  25 
November  £uid  29  September  1987.  The  data  correspond  to  low  and  moderate 
current  wind  speeds  )  of  2.5  m/s  and  7.5  m/s  for  the  respective  days.  The 
Thermal  Video  Processor  System  (THERMOTEKNIX)  available  with  the  AGA 
system  also  allows  the  thermal  scene  to  be  displayed  in  a  format  consisting  of  128 
pixel  lines  (0.023  degAine).  The  effective  blackbody  temperature  corresponding 
to  each  pixel  can  then  be  displayed  on  the  screen  by  positioning  a  cursor  at  the 
appropriate  position.  Vertical  temperature  profiles  in  the  scene  are  displayed  on 
the  right-side  of  the  thermogram  as  shown  in  Fig.  la  and  2a.  In  each  case,  the 
horizontal  cursor  is  situated  on  the  pixel  corresponding  to  the  maximum  tem¬ 
perature  (20.4 °C  or  3.46  mW/cm^  sr  and  16.5®C  or  3.23  mW/cm*  sr)  which  is 
taken  to  coincide  with  the  infrared  horizon.  The  increase  in  temperature  with 
altitude  (decreasing  zenith  angle)  in  the  29  September  scene  reflects  the  increeise 
in  air  temperature  with  altitude  shown  in  the  Fig.  3.  Similarly,  the  decrease  in 
temperature  with  altitude  in  the  25  November  scene  can  be  related  to  the  de¬ 
crease  in  air  temperature  with  altitude.  The  larger  sky  radiance  values  on  29 
September  may  also  partiedly  result  from  the  lower  relative  humidity  above  200 
m  (Fig.  4),  which  decreeised  the  number  of  activated  aerosols  available  to  scatter 
the  radiation. 

The  THERMOTEKNIX  system  also  allows  the  relative  frequency  of  occur¬ 
rence  of  pixels  with  a  specific  temperature  to  be  calculated  for  selected  rectangu¬ 
lar  areas  within  a  scene.  When  the  histogram  is  plotted,  a  Gaussian  curve 
(shown  with  artist-enhanced  dots)  is  calculated  which  best  fits  the  temperature 
points.  In  both  the  low  and  moderate  wind-speed  examples,  the  sIq^  temperature 
distributions  (Fig.  lb  and  2b)  are  more  closely  Gaussian  than  are  those  of  the  sea 
(Fig.  Ic  and  2c)  which  show  a  banded  feature.  The  mean  temperature  differences 
for  the  two  examples  are  3.6®C  and  1.4°C,  respectively.  In  the  higher  wind  speed 
(9.3  m/s)  example  shown  in  Fig.  5  for  19  March  1987,  when  whitecaps  were  pre¬ 
sent,  the  sea  temperature  distribution  (Fig.  5c)  is  remarkably  similar  to  a  Gauss¬ 
ian  distribution.  Tbe  mean  temperature  difference  between  the  sky  and  sea  in 
this  case  is  only  0.2 °C.  In  contrast,  the  sea  and  sl^  mean  temperature  distribu¬ 
tions  shown  in  Fig.  6  during  stratus  cloud  conditions  (29  April  1987)  were  found 
to  be  equsd,  euid  both  were  closely  Gaussieui.  The  surface  wind  speed  was  not 
measured  in  this  example. 
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For  this  study,  18  thermograms  recorded  during  different  meteorological 
and  wind-speed  conditions  in  1987  and  1988  were  used  to  compare  the  mean 
cloud-free  sky  and  sea  temperatures.  In  Table  1,  the  measured  mean  sky  and  sea 
temperatures  and  the  current  wind  speed  are  shown.  In  12  of  the  cases,  the 
24-hour  averaged  wind  speed  (  F)  and  the  actual  sea  surface  temperature  (r,,) 
measured  by  the  aircraft  were  available.  Also  shown  are  the  combinations  of  air 
mstss  factors  (AM)  (see  Appendix  B)  and  surface  visibilities  (VIS),  which  allow  the 
horizon  pixel  radiance,  c^culated  using  LOWTRAN  6  (with  the  measured  pro¬ 
files  of  meteorological  parameters)  to  exactly  match  the  measured  value  (Hughes 
euid  Jensen,  1988).  For  this  data  set,  the  mean  temperatures  were  highly  corre¬ 
lated  (correlation  coefficient,  r  =  0.92),  as  shown  in  Fig.  7,  with  the  sea  tempera¬ 
ture  being  less  than  the  sky  temperature,  as  indicated  by  the  dashed  line  for 
one-to-one  correspondence.  The  mean  temperatures  differed  the  most  during  low 
wind-speed  conditions.  As  shown  in  Fig.  8,  the  differences  were  found  to  de¬ 
crease  with  increasing  wind  speed. 

CALCULATION  OF  BACKGROUND  SCENES 

Earlier  work  (Hughes  and  Jensen,  1988)  has  demonstrated  that  the  sky 
radiances  measured  at  individual  pixel  lines  within  about  1  degree  of  the  horizon 
can  be  modeled  very  closely  by  using  the  LOWTRAN  6  code.  Figures  9  and  10 
show  the  comparison  of  the  measured  and  calculated  radiances  for  zenith  angles 
within  about  1  degree  below  the  horizon  for  the  scenes  in  Fig.  1  and  2.  In  the 
calculations,  the  meteorolopcal  profiles  were  divided  into  33  layers,  as  allowed 
by  LOWTRAN  6.  The  lower  layers  of  the  profiles  were  also  divided  into  sub¬ 
layers  containing  the  same  amount  of  absorbing  and  scattering  material  and  the 
temperature  of  the  original  layer.  This  artificial  layering  has  been  found  neces¬ 
sary  (Wollenweber,  1988a)  to  remove  the  anomolous  dip  (Hughes,  1987)  which 
occurs  when  aerosols  are  included  in  the  radiance  calculations  for  zenith  angles 
close  to  90  degrees.  As  the  AGA  scanner  could  not  be  accurately  plumbed,  the 
zenith  angle  of  the  infreu'ed  horizon  (the  pixel  corresponding  to  the  maximum  ra¬ 
diance)  in  each  thermogram  was  teiken  to  be  one-half  of  a  pixel  less  than  the  an¬ 
gle  (0  =  90.179  degrees  on  29  September  and  0  =  90.174  degrees  on  25  Novem¬ 
ber)  at  which  the  LOWTRAN  calculations  indicated  the  refracted  ray  path  first 
struck  the  e6ui.h.  In  both  cases,  the  meyor  contributor  to  the  total  radiance  just 
below  the  horizon  is  the  path  emission  N{e)p  .  While  the  reflected  sky  radiance 
N(0)„k  and  the  surface  emission  N{0)ss  are  small,  their  contribution  to  the  total 
radiance  N{G),ot  at  this  low  level  (33  m)  of  observation  cannot  be  neglected.  It  is 
interesting  to  note  only  a  small  reversal  of  the  relative  magnitudes  of  the  re¬ 
flected  sl^  radiances  and  surface  emissions  between  the  two  sets  of  calculations, 
which  demonstrates  the  small  influence  of  the  wave  slopes  for  the  moderate  wind 
speeds  on  29  September. 

In  a  similar  manner,  the  mean  equivalent  blackbody  temperatures  corre¬ 
sponding  to  1-degree  eu'eas  above  and  below  the  horizon  were  calculated  for  the 
12  cases  where  vertical  profiles  of  meterological  peu’ameters  were  available.  In 
Fig.  11,  the  meeisured  and  calculated  mean  sky  temperatures,  Tm(,Sky)  and  T^Sky), 
are  in  close  eigreement  (r  =  0.99).  In  Fig.  12,  the  calculated  mean  sea  tempera¬ 
tures,  Tc(Sea),  show  slightly  more  variations  with  the  measured  values,  T„(Sea) , 
but  the  two  are  still  fairly  well  correlated  (r  =  0.90). 
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DISCUSSION 


The  demonstrated  agreement  between  the  measured  and  calculated  mean 
temperatures  places  confidence  in  the  usefulness  of  the  models  to  calculate  radi¬ 
ances  as  viewed  from  other  altitudes  In  a  practical  sense,  the  sea  radiance 
model  needs  to  be  improved  to  reduce  computer  running  time  In  its  present 
form,  approximately  3  minutes  of  computational  time  on  the  HP-9020  are  re¬ 
quired  to  calculate  the  sea  radiance  at  one  zenith  angle  of  viewing  as  compared 
to  approximately  5  seconds  for  the  sky  radiance  In  this  study,  approximately  30 
minutes  of  computer  time  was  required  to  calculate  the  mean  sea  temperature 
within  I  degree  of  the  horizon  on  each  day.  The  m^or  contributor  to  the  compu¬ 
tational  time  is  the  myriad  of  sky-radiance  reflections  and  surface-emission  cal¬ 
culations  that  must  be  made  to  account  for  the  radiance  from  all  the  wave  facets 
C/omputational  time  could  be  saved  if  the  surface  emissivity  and  reflectivity  cal¬ 
culations  were  made  at  a  single  representative  wavelength  instead  of  having  to 
be  averaged  over  a  wavelength  band  for  each  wave  slope  Another  way  of  reduc¬ 
ing  the  computational  time  is  to  develop  an  empirical  relationship  to  directly  re¬ 
late  the  apparent  mean  sea  temperature  to  the  calculated  mean  sky  temperature. 

Inherent  in  the  observed  variations  between  the  mean  sky  and  sea  tem¬ 
perature  differences  in  Fig.  5  are  the  surface  wind  speeds  and  actual  sea  surface 
temperatures.  If  these  pareimeters  are  included  in  a  multiple  regression  analysis 
for  three  independent  variables,  we  obtain  the  following  relationship: 

r«(5ca)  =  1.09  T„(Sky)  +  0.37  ^  0.24  T„  -  10.46  (1) 

The  correlation  coefficient  for  the  multiple  regression  analysis  is  0.97.  Only  12 
samples  were  used  in  the  analysis  (as  compared  to  18  samples  in  the  linear  re¬ 
gression  analysis  shown  in  Fig.  5).  which  indicates  a  definite  improvement  in  the 
correlation  when  the  actual  sea  surface  temperatures  and  wind  speeds  are  in¬ 
cluded. 


Based  on  earlier  work  (Hughes  and  Jensen,  1988),  which  demonstrated 
that  the  sky  radiances  measured  at  individi  al  pixel  lines  within  about  1  degree  of 
the  horizon  can  be  modeled  very  closely  by  using  the  LOWTRAN  6  code,  and  the 
excellent  agreement  between  the  calculated  and  measured  mean  sIq^  radiances  in 
Fig.  11,  the  replacement  of  T„iSky)  by  a  calculated  value  of  the  mean  sky  tem¬ 
perature,  Tc(.Sky) ,  would  appear  justified.  This  would  allow  the  mean  sea  tem¬ 
perature  to  be  directly  determined  by  using  Eq.  1  and  measured  profiles  of 
meteorological  pcu-ameters.  Additional  data  are  needed  to  verify  that  such  a  rela¬ 
tionship  for  predicting  the  apparent  mean  temperature  of  the  sea  is  meaningful 
for  other  meteorological  conditions,  sea  states,  and  altitudes  of  observation. 

Such  a  relationship  would  be  a  valuable  inclusion  into  a  system  designer’s  hand¬ 
book  as  a  tool  for  predicting  infrared  sea  backgrounds  in  locales  where  only  mete¬ 
orological  data  are  available. 
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Appendix  A 

MATHEMATICAL  FORMULATION  OF  SEA  RADIANCE 

MODEL 

Consider  the  atmosphere  to  be  composed  of  a  number,  n,  of  isothermal 
layers  characterized  by  temperature  Ti  and  transmittance  T(v,i,n)  along  the  opti¬ 
cal  path  traversing  the  tth  layer  at  an^c  n,  and  u  is  the  spectral  wave  number. 
From  KirchofTs  law,  the  radiance  of  the  ith  layer  is 

N(u,j,n)j*  =  [1  -  r<j  (ii,i,n)]W(r,  Vtt  (A-1) 


where  (w,i,ji)  is  the  absorption  transmittance  and  W(T,  )  is  Planck’s  blackbody 
radiation  formula.  Then  the  spectral  radiance  reaching  the  sea  surface  throu^ 
the  intervening  atmosphere  is 


N(v,i,ii)sk  n  r(v,;./i)  =  (1  -  rfl(u,i,n)]  U  T(v,j,fi)  W{Ti)h 

_y  =  1  J  > 


(A-2) 


Summing  the  contribution  from  all  layers,  the  spectral  radiance  at  the  sea  sur¬ 
face  is  then 


N{v,\y)sk  = 


Mv,i,p)]r'n  r(v.y.^)  W{Ti)h) 

)  =  i  U*> 


(A-3) 


The  radiance  is  allowed  to  strike  a  wave  facet  on  the  ocean  surface  with  a  Gauss¬ 
ian  distribution  (Cox  and  Munk,  1954)  of  angular  tilts  a  and  p  in  the  up-wind 
and  cross-wind  directions,  respectively,  so  that  £Ui  amount  N(v,>i)'  is  reflected 
into  the  sensor  at  an  edtitude  //,  within  the  mth  layer.  The  probability  that  radi¬ 
ance  hits  the  facet  is  equ£d  to  the  probability  that  the  wave  slope  exists,  i.e., 

N(v,n)'  =  PiS,,Sy)N{v,v.)sk  (A-4) 


where 


P{Sx,Sy)  =  l/(2a;,ay)EXP[0.05(5iM  (A-5) 

and  Sx  =  tana,  Sy  =  tanp,  oi  =  0.003  +  1.92xl03V^c,  oy  =  3.16xl03Vc,  with  Vc 
being  equal  to  the  current  wind  speed  in  the  azimuthal  direction  4)  with  respect 
to  the  sensor.  Then  the  total  spectral  radiance  that  is  reflected  from  all  the 
wave  facets  into  the  line-of-sight  of  the  detector  located  in  the  mth  layer  is 


Niv,Q)rsk 


"'u  T(j,  e) 


;  =  > 


2  R 


iv,n)P{Sx  ,  Sy)N{v,ik)rsk 


(A-6) 


where  is  the  complex  reflectivity  of  seawater  at  the  reflection  an^e  n.  In 

the  above  equations,  both  n  and  n  are  implicit  functions  of  Sx  and  Sy  given  by 
(Wollenweber,  1988b) 


5 


cos^  =  (25;t/A)cos0'cos4>  +  (2Sy/A)cos0'cos<|>  -  (B/A)sin0' 
cosn  =  (5x/A)cos0'cos<t))  +  (i>/A)cos0'cos<j>  +  (1/A)sin0' 


(A-7) 


(A-8) 

where  A  =  Sj  +  +  1,  B  =  Sl  +  Sy  -  1,  and  0'  is  the  sensors’  zenith  6Ln|^e  at 

the  sea  surface  reflection  point. 

In  a  similar  manner,  the  spectral  radiances  emitted  by  the  sea  surface 
wave  facets  (Nss)  ar  1  the  path  radiance  which  reach  the  sensor  at  the  zenith 
angle  0  are  given  by 

N{v,e)ss  =  r(v.y.  0)  j  ^(5;,  ,  Sy)[l  .  /Jfv.nlimT-wl/TT  (A-9) 


and 


N(v,Q)p  =  X  [1  -  (i’,i,0)] 


I  =  1 


m-1  I 

n  T{vj,e) 

y  = '  ♦  1  J 


WiTi)h 


(A-10) 


where  Tss  is  the  sea  surface  temperature,  and  again,  the  an^e  g  is  implicit  in  the 
reflection  angle  n. 

Then,  the  total  spectral  radiance  reaching  the  detector  is  the  sum  of  the 
three  components: 

N{v,9)toi  -  Niv,Q)rsk  +  N(v,Q)ss  +  N{v,9)p  (A-11) 

The  total  spectral  radiance  must  then  be  averaged  over  the  response  of  the  sensor 
system,  which  in  this  case  is  taken  to  be  the  8-  to  12-fim  wavelength  band.  Sub¬ 
routines  have  been  introduced  into  LOWTRAN  6  (Wollenweber,  1988b)  to  calcu¬ 
late  the  band-averaged  radiance.  The  reflection  and  zenith  an^es  are  calculated 
with  Eq.  A-5,  A-7,  and  A-8,  corresponding  to  the  incremented  values  of  wave 
slopes  in  the  intervals  -Zox.y  <  Sx,y  <  Sox.y. 
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Appendix  B 

DESCRIPTION  OF  THE  LOWTRAN  6  NAVY  MARITIME 

AEROSOL  MODEL 

The  peu-ticle  size  distribution  model  (at  radius  r 
normal  distributions  given  by 

)  is  the  sum  of  three  log- 

3 

n{r)  =  X  ^  1  exp[-fln )2]  (cm-3  •  pjn-l) 

(B-1) 

where 

A,  =  2000(AM)2 

(B-2) 

A2  =  5.866(F  -  2.2) 

(B-3) 

A3  =  10  {OMVc  -  2.8) 

(B^) 

Component  A,  represents  the  contribution  by  continental  aerosols.  AM  is 
an  air  mass  parameter  that  is  allowed  to  range  between  integer  values  of  1  for 
open  ocean  and  10  for  coastal  areas  and  is  given  by 

AM  =  INTiRnIA)  +  1  (B-5) 

where  Rn  is  the  measured  atmospheric  radon  content  expressed  in  pCi/m^.  In 
the  absence  of  radon  measurements,  the  air  mass  factor  can  be  related  to  the 
elapsed  time,  T(days),  for  the  air  mass  to  reach  the  point  of  observation: 

AM  =  /^T[9exp(-T/4)]  +  1  (B-6) 

Components  A2  and  Aa  represent  equilibrium  sea  spray  particles  generated 
by  the  surface  wind  speed  averaged  over  24  hours  (V,  in  m/s)  and  the  current  sur¬ 
face  wind  speed  (Vc,  in  m/s),  respectively.  (It  should  be  noted  that  the  current 
wind  speed  component  is  different  from  the  value  published  in  LOWTRAN  6. 

This  modification  was  found  to  be  necessary  in  order  to  match  previously  pub¬ 
lished  measurements  of  infrared  slqr  radiances  and  near-surface  aerosol  size  dis¬ 
tributions  (Hughes,  1987)).  In  Eq.  B-1,  r,- ,  the  modal  radius  for  each  component 
referenced  to  a  relative  humidity  of  80%  (r,  =  0.03  >i.m,  Tz  =  0.24  >im,  and  = 
2.0  Jim),  is  allowed  to  grow  with  relative  humidity  {RH)  according  to  the  formula 
(Fitzgerald,  1975) 

/•  =  r(2  -  /?H/100)/6(1  -  i?if/100)lV3  (B-7) 

The  contribution  to  the  total  extinction  or  absorption  by  each  component 
can  then  be  written  u*. 


where  Ci  =  (O.OOIW/)-^/. 
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The  factor  in  the  expression  for  Ci  ensures  a  constant  total  number  of 
particles  as  the  relative  humidity  increases.  Qt.a  is  the  cross-section  for 

either  the  extinction  or  absorption  normalized  to  the  geometrical  cross-section  of 
the  spherical  particle,  and  m  is  the  complex  refractive  index,  which  is  allowed  to 
change  from  that  of  dry  sea  salt  as  the  particle  deliquesces  with  increasing  hu¬ 
midity.  LOWTRAN  6  provides  precalculated  values  in  tabular  form  of  the  pa¬ 
rameter  ae,a  {k)l Ci  at  discrete  wavelengths  for  four  relative  humidities  (50%, 
85%,  90%,  and  99%),  from  which  the  average  extinction  for  a  specific  wavelength 
band  and  relative  humidity  can  be  readily  determined  by  interpolation.  When  an 
observed  surface  visibility  {VISo)  is  available  as  an  input  to  the  model,  the  ampli¬ 
tudes  of  the  three  components  are  adjusted  by  a  scaling  factor  (SF)  so  that  the 
calculated  aerosol  extinction  coefficient,  at  a  wavelength  of  0.55  pm  is  the 
same  as  the  observed  extinction,  Oo,  determined  from  the  relationship 


VISo 


3.912 

ao*ar 


(B-9) 


where  Or  is  the  Rayleigh  contribution  to  extinction  at  0.55mm. 
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LIST  OF  CAPTIONS 


Figure  1 .  'mermograms  of  near-horizon  infrared  radiances  on  25  November  1987  during  low 
wind  speed  conditions  (2.S  nVs).  (a)  Vertical  profile  of  effective  biackbody  tempera¬ 
ture  near  the  horizon,  and  temperature  histograms  coaesporxling  to  areas  1  de^ee 
above  (b)  and  1  degree  below  (c)  the  horizon,  m  each  histogram  a  Gaussian  curve 
is  displayed  as  a  series  of  dots  which  best  fit  the  temperature  points. 

Figure  2.  Thermograms  of  near-horizon  infrared  radiances  on  29  September  1987  during  mod¬ 
erate  wind  speed  conditions  (7.5  m/s)  without  the  presence  of  whitecaps.  (a)  Verti¬ 
cal  profile  of  effective  biackbody  temperature  near  the  horizon,  and  temperature 
histograms  corresponding  to  areas  1  degree  above  (b)  and  1  degree  below  (c)  the 
horizon.  In  each  histogram  a  Gaussian  curve  is  displayed  as  a  series  of  dots  which 
best  fit  the  temperature  points. 

Figure  3.  Profiles  of  air  temperature  measured  with  altitude  on  29  September  and  25  Novem¬ 
ber  1987  off  the  coast  of  San  Diego,  California. 

Figure  4.  Profiles  of  relative  humidity  measured  v4th  altitude  on  29  September  and  25  Novem¬ 
ber  1987  off  the  coast  of  San  Diego,  California. 

Figure  5.  Thermograms  of  near-horizon  infrared  radiances  on  19  March  1987  during  moderate 
wind  speed  conditions  (9.2  m/s)  when  whitecaps  were  present,  (a)  Vertical  profile  of 
effect!^  biackbody  temperature  near  the  horizon,  arxl  temperature  histograms  corre¬ 
sponding  to  areas  1  degree  above  (b)  and  1  degree  below  (c)  the  horizon.  In  each 
histogram  a  Gaussian  cunre  is  displayed  as  a  series  of  dots  which  best  fit  the  tem¬ 
perature  points. 

Figure  6.  Thermograms  of  near-horizon  infrared  radiances  on  29  April  1987  during  stratus  cloud 
condHions  (wind  speed  not  measued).  (a)  Vertical  profile  of  effective  biackbody 
temperature  near  the  horizon,  and  temperature  histograms  corresporKling  to  areas  1 
degree  above  (b)  and  1  degree  below  (c)  the  horizon.  In  each  histogram  a  Gaussian 
cunre  is  displayed  as  a  series  of  dots  which  best  fit  the  temperature  points. 

Figure  7.  Regression  analysis  of  measured  mean  sky  temperatures,  7m(5ity)>  ahd  mean  sea 
temperatures,  Jm  {Sea}  The  dashed  line  corresponds  to  a  one-to-one  correspon¬ 
dence  between  the  two  parameters. 

Figure  8.  Regression  analysis  of  the  difference  between  the  measued  mean  sky  temperatues, 
TmiSky).  and  mean  sea  temperatues,  TmiSea}  and  the  surface  wind  speed. 

Figue  9.  Comparison  of  the  measued  and  calculated  IR  radiances  fOr  zenith  angles  about  1 
degree  below  the  horizon  on  25  November  1987. 

Figue  1 0.  Comparison  of  the  measued  arxl  calculated  IR  radiances  for  zenith  angles  about  1 
degree  below  the  horizon  on  29  September  1987. 

Figue  1 1 .  Coaelation  between  the  measued,  T„(Siy}  and  calculated,  Te(S*y) .  niean  sky 
temperatues. 

Figue  12.  Correlation  between  the  measued,  T„(Sea) .  and  calculated,  Tc(Sea) .  mean  sea 
temperatues. 
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TABLE  1.  Mean  sky,  TmiSky) ,  and  sea,  Tm(Sea\  temperatures,  current  wind  speed,  Vc, 

24-hour  averaged  wind  speed,  and  the  actual  sea  surface  temperature,  T,s,  measured  on 
the  days  indicated.  VIS  and  AM  are  the  inferred  visibility  and  air  mass  factor,  respectively. 


Date 

T,(Sky) 

(°C) 

T„(Sea) 

(°C) 

Vc 

(m/s) 

V 

(m/s) 

^ss 

(°C) 

VIS 

(km) 

(1987) 

19  Mar 

10,1 

9.9 

9.3 

25  Mar 

9.5 

4.6 

1.5 

26  Mar 

12.9 

8.9 

4.1 

3  Apr 

11.2 

10.5 

8.8 

18  May 

14.3 

12.8 

6.1 

26  May 

13.6 

13.3 

11.6 

30  Jun 

16.4 

14.3 

68 

5.0 

17.8 

45.0 

8  Sep 

21.6 

18.8 

4  7 

2.6 

18.7 

16.5 

17  Sep 

19.0 

16.9 

69 

2.9 

18.1 

31.5 

29  Sep 

20.3 

18.9 

7.5- 

4.4 

15.2 

26.5 

3  Nov 

16.0 

13.7 

5.3 

3.9 

17.7 

24.6 

12  Nov 

18.6 

15.4 

3.2 

2.7 

17.0 

24.5 

24  Nov 

16.7 

13.8 

[ 

6.8 

3.4 

17.4 

18.0 

25  Nov 

16.1 

12.5 

2.5 

3.9 

17,6 

21.0 

(1988) 

9  Jun 

16.1 

1 

12,1 

2.9 

2.8 

16.4 

37.0 

10  Jun 

16.1 

12.1 

4.6 

2.2 

16.8 

35.5 

9  Nov 

16.7 

13.6 

3.7 

2.1 

19.3 

41.3 

10  Nov 

17.0 

15.6 

4.8 

2.8 

20.3 

45.8 

